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dimensionless potential temperature gradient 
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TECHNICAL MEMORANDUM X- 64519 


THE POTENTIAL TEMPERATURE PROFILE IN THE PLANETARY BOUNDARY LAYER 

SUMMARY 


The profile of mean flow potential temperature in the first 
150 meters of unstable planetary boundary layer is analyzed with a 
sample of nineteen temperature profiles observed at the NASA 150-meter 
meteorological tower site. The observations consist of mean tempera- 
tures observed at 3-, 18- , 30- , 60- , 120- , and . 150-meter levels. The 
duration time of each test ranged between approximately fifteen minutes 
to one hour. 

It is hypothesized that the dimensionless potential temperature 
gradient 0q = zT^o^B/^z is a dimensionless function of z, L Q , f, and 
u* 0 , where §(z) is the mean flow potential temperature at height z, 

T* 0 and u* 0 denote the surface values of the friction temperature and 
friction velocity, L 0 is, the surface Monin-Obukhov stability length, 
and f is the Coriolis parameter. According to Buckingham's theorem, 
only two independent dimensionless quantities can be constructed from 
the set z, L 0 , f and u* Q ; however, these dimensionless quantities can 
be chosen in a variety of ways. Accordingly, the dependencies of 0q 
on (z/L Q , R-^)-, (z/L c , fz/u* Q )-, and (fz/u* 0j R^) -coordinates are 
examined, where Rj, = -fL 0 /u* 0 . The function 0q(z/L q , R l ) = (l-yz/L 0 )“ x 
summarized the data reasonably well, where 7 is a positive quantity pro 
portional to R£ 3 / 2 over the range of variation of the data. The func- 
tions 0q( z/L q , fz/u* Q ) and 0Q(fz/ti* o , R^) , derived from the function 
^0(z/L o , R^) , agree with the experimental results. 

Integration of the differential equation 0 g = 0 q(z/L o ,R l ) yields 
the potential temperature profile. For* sufficiently small values of 
the quantity A = -7 z 0 /L q , the potential temperature profile has a 
logarithmic behavior with a linear defect resulting from stability and 
Coriolis effects, where z Q is the surface roughness length of the site. 
The potential temperature profile is used to evaluate the drop in 
potential temperature across the planetary boundary layer. The depth 
of the boundary layer is assumed to be u >vo /4f. 



I. INTRODUCTION 


This report concerns the vertical profile of mean flow potential 
temperature in the planetary boundary layer. In approximately the first 
30 meters of the horizontally homogeneous planetary boundary layer, the 
Monin layer, the experimental results of Dyer [1], Pandolfo [2], and 
Prasad and Panofsky [3] imply that the potential temperature profile can 
be expressed functionally in the form 


0 F (z/L q , z o/Lq) > 


( 1 ) 


where z, L 0 , and z Q denote the height, surface Monin-Obukhov stability 
length, and the surface roughness length. The mean flow wind profile is 
also a function of z/L 0 and z 0 /L 0 in the Monin layer. However, as we 
proceed upward out of the Monin layer into the Ekman layer, we find that 
the Coriolis forces must be taken into account in order to explain the 
behavior of the wind profile [4, 5]. In view of the inherent coupling 
between the vertical turbulent momentum and heat fluxes and the mean 
flow wind and potential temperature profiles implied by the mean-flow 
Boussinesq-approximated equations of motion [6], it is reasonable to 
expect that the potential temperature profile is also affected by the 
Coriolis terms in the momentum conservation equations. The data analysis 
that follows shows that this appears to be the case. 


II . DIMENSIONAL CONS IDERATIONS 


In the Monin layer the vertical gradient of mean flow potential tem- 
perature 50/Bz is scaled with the height z and the surface friction tem- 
perature T* 0 . The resulting dimensionless temperature gradient is a 
universal function of z/L 0 , so that 


t* q dz 0 9 (z/l o ). 


( 2 ) 


The quantity L 0 is the surface Monin-Obukhov stability length and is given 
by ' 


C P P T u | 0 

V ’ k lg H q ’ 


( 3 ) 
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where u^ Q is the surface friction velocity, k ± is von Karman f s constant 
with a numerical value approximately equal to 0,4, T and p denote the 
mean temperature and density in the boundary layer, g is the accelera- 
tion of gravity, and H Q is the surface heat flux. The surface friction 
temperature is defined as 


H q 

c p ' 


(4) 


Integration of equation (2) from z Q to z, subject to the boundary condi- 
tion that 


0 = 0(z o /L o ) at z/L q » z 0 /L 0 , 


(5) 


will yield the potential temperature profile 


0(z/L o ) - 0(z o /L o ) + T* 0 



1 - 0 e (x) 

X 



( 6 ) 


Note that the surface roughness length appears in the potential tempera- 
ture profile through the boundary condition (5), 

In the Ekman layer, 0q should also be a function of z and L 0 and 
other parameters that characterize the action of Coriolis forces, 
baroclinic effects, etc. In this report, we are concerned with deter- 
mining the effects of the vertical heat and momentum fluxes and the 
Coriolis forces on the potential temperature profile. The parameter 
L 0 should be sufficient to characterize the effects resulting from the 
turbulent momentum and heat fluxes. To represent the effects of 
Coriolis forces, we could add the Coriolis parameter f to our list of 
independent variables,* Thus, we might suspect that 0q should depend 
on z, L 0 , and f; however, it is not possible to construct a dimension- 
less quantity that contains f with this list of independent variables, 
because z and L 0 have the units of length and f is an inverse time. 


The Coriolis parameter is defined as f = 2q sin 0 , where Q is the 
angular velocity of the earth and 0 denotes the latitude of the site. 
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Thus, at least one additional parameter must be included in the dimen- 
sional analysis, Blackadar and Tennekes [4] have shown that u* Q is the 
appropriate velocity scale in the planetary boundary layer by analyzing 
the turbulent energy equation. The addition of u^q to our list of 
pertinent parameters will enable us to construct a dimensionless 
quantity which will contain f. This follows from the fact that u* Q / f 
has the units of length. In fact, an estimate of the thickness of the 
planetary boundary layer is u..yo/4-f [4], so that the addition of u*o to 
our list of independent variables is physically relevant. As in the 
case of the Monin layer, we will not include z Q in the list of variables 
upon which we hypothesize 0§ to depend. However, z Q will enter the 
analysis through the boundary condition, equation (5), when we integrate 
0 q to obtain the potential temperature profile. 

According to Buckingham’s theorem [7], the number of a complete set 
of independent dimensionless quantities that can be constructed from the 
set of four variables z, L 0 , f, and u v?0 is two. There are many possible 
ways to construct these dimensionless quantities. Three possible repre- 
sentations are given by 


Set I: 


Set II: 


Set III: 


z 

fL ■> 
o ! 


L 5 
o 

u * 0 ! 


z 

f z 


L 5 
o 

u *o 

• 

f z 

fL 
0 


u *o’ 

u *o j 



(7) 


These dimensionless quantities are the simplest representations since 
they are constructed from the independent variables raised only to the 
first power. All other formulations would require lower or higher order 
exponents. Any one set of the dimensionless quantities given by (7) can 
be derived from ahy one of the remaining two sets by forming appropriate 
ratios. Each set has its own merits so that, rather than selecting one 
formulation of the problem, we will examine the following three repre- 
sentations of 0 q : 


0 Q = 0 e ( Z /L o , R l ) 


( 8 ) 
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0e 


0 e (z/L o , fz/u* Q ) 


(9) 


where 


00 = 0 e (fz/u* o . r l>» 


fL, 


*L = ' u 


*o 


( 10 ) 


(ID 


III. THE DATA SOURCE AND DATA PROCESSING 


The data analyzed in this report consist of nineteen sets of tem- 
perature measurements at the 3-, 18- , 30- , 60- , 120- , and 150-meter 
levels and associated wind speed measurements at the 18- and 30-meter 
levels obtained at NASA's 150-meter meteorological tower site. The 
instrumentation of this tower site has been discussed by Kaufman and 
Keene [8], The azimuthal dis tribution .of the surface roughness z Q at 
the KSC tower, site is discussed in reference 9, and the temporal mean 
temperature and wind data are' tabulated in Appendix B. The duration of 
each case ranged between approximately fifteen minutes to one hour; how- 
ever, the majority of cases had a duration time of one hour. The cal- 
culation of u* 0 , T* 0 , and L 0 for each case is discussed in Appendix A, 
and the results of the calculations are tabulated in Appendix B. 

The mean flow potential temperature 0, Kelvin temperature T, and 
pressure p are related through Poisson's equation 


0 = T^/p) 



( 12 ) 


where Pi = 1,000 mb and R is the specific gas constant for dry air. 
Partial differentiation of this equation with respect to z yields the 
result 


SZ Ti 



(13) 


Equation (13) was used to obtain a first estimate of Si/Bz, with the 
assumption that 0/T ~ 1. The derivative ST/^z was estimated at the 
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mid-points between the data acquisition levels with finite central dif- 
ferences. The resulting estimates of Bo/Bz were multiplied by correc- 
tion factors to account for the errors resulting from the finite 
difference operators. These corrections were based on the assumption 
that Be/Bz « z^ over piecewise portions of the potential temperature 
profile, where q is a constant. The details of this correction and 
its application are discussed in Appendix C. 


IV. DIMENSIONLESS POTENTIAL TEMPERATURE GRADIENTS 


The ranges of the experimental values of R^ and fz/u, vo were 
0.0004 < Rl < 0.05 and 0.0007 < fz/u*o < 0.03, The nineteen poten- 
tial temperature profiles were grouped according to the categories of 
R^ given in Table I and the ninety-five estimates of the dimensionless 
potential temperature gradient were grouped according to the categories 
of fz/u* Q given in Table II. Figures 1 through 3 are plots of the 
experimental values of 0 0 (z/L o , R L ), 0 0 (z/L o , fz/u* 0 ), and 0 0 (f z/u* 0> Rl) 
with the data points grouped according to the categories in Tables I or II. 


TABLE I 

Potential Temperature Profile Categories According to R^ 
Category No. of Profiles 


0. 

,0004 

< 

r l 

< 

0. 

,0007 

0 . 

,001 

< 

*L 

< 

0. 

,003 

0 . 

,005 

< 

r l 

< 

0 . 

,01 

0. 

,01 

< 

*L 

< 

0. 

,03 

0. 

,04 

< 

r l 

< 

0 . 

,05 


TABLE II 

Potential Temperature Profile Categories According to fz/u, v0 

No. of Qbs er va t ions 


0 . 

,0007 

< 

fz/u*o 

< 

0. 

,001 

4 

0 . 

,001 

< 

fz/u*o 

< 

0. 

,002 

16 

0. 

,002 

< 

fz / u *o 

< 

0 . 

,004 

17 

0 . 

,004 

< 

fz/u*o 

< 

0 . 

,006 

14 

0. 

,006 

< 

fz/u*o 

< 

0 . 

,01 

15 

0 . 

,01 

< 

fz/u*o 

< 

0 . 

,03 

29 
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Figure 1 shows that 0q is an increasing function of z/L 0 for fixed 
Rl and an increasing function of Rj, for fixed z/L q . It might be con- 
cluded from these results that the presence of Coriolis forces tends to 
promote stability in the sense that an increase in f will tend to pro- 
duce an increase in S 0 /Sz, all other things being equal. 

Figure 2 shows that 0q is an increasing function of z/L 0 for fixed 
fz/u * 0 and an increasing function of fz/u * 0 for fixed z/L 0 . Figure 3 
implies that 0q is a decreasing function of fz/u * 0 for fixed R^ and an 
increasing function of R^ for fixed fz/u, vo . The upper bound value on 
00 is of order one, while the lower bound appears to approach zero like 
(z/L 0 ) -1 / 2 for fixed R-r in figure 1, (z/L 0 )" 5 / 4 for fixed fz/u* D in fig*“ 

ure 2, and (f z/u * 0 ) -1 ' 2 for fixed R^ in figure 3. 

The Businger hypothesis states that in the Monin layer the flux 
Richardson number 


& Se 

Ri - - Q — — (14) 

(Su/Sz ) 2 


is related to L Q through the expression 


Ri = z/L OJ (15) 

where u is the mean wind at height z. Experimental investigations [ 6 ] 
imply that the dimensionless wind shear 0 u = (k 1 z/u, v o)Su/Sz in the Monin 
layer is given by 


0 u = (1 - 18Ri) -1 / 4 . 


(16) 


Upon combining (3), (4), (14), 


(15), and (16), we find that 


0 Q = (1 - 18z/L 0 ) -1 / 2 . 


(17) 


Thus, as “z/L 0 becomes large in the unstable Monin layer, 0q asymptoti- 
cally behaves like 


00 ~ 18 -1 / 2 (z/L 0 ) -1 / 2 . 


(18) 
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We should be able to obtain this asymptotic behavior from our data as 
Rl -> 0. The experimental values of 0 q in Figure 2 tend to decrease as 
(-z/Lo )* 1 / 2 as -z/L 0 becomes large for the smaller values of How- 

ever, the coefficient of (-z/L 0 )“ 1 / 2 in equation (18) is not obtained 
as Rl becomes small, but rather the data appear to imply that the coef- 
ficient increases to values greater than 18" 1 / 2 as R L becomes small* 
This will become more obvious below. 

To develop a model of 0 q for the planetary boundary layer, it is 
assumed that 

J2f 0 (z/L O ,R L ) = (1 - /z/Lo)- 1 / 2 , (19) 


where 7 is now a function of Rl* When 7 = 18, we obtain equation (17). 
The function 7 was determined by fitting equation (19) to each potential 
temperature profile with least-square procedures. The results of these 
computations are shown in Figure 4. Three of the profiles yielded nega- 
tive values of 7 which were approximately equal to zero for R^ = 0.0237, 
0.0415, and 0.0474. However, these data points fit in with the overall 
trend implied- by the data points in Figure 4, namely, 7 being a decreas- 
ing function of R L . A least-squares analysis of the sixteen data points 
in Figure 4 yielded the result 


7 = 0.0044 R£ 3 ^ 2 . 


( 20 ) 


It is of interest to note that the geometric mean value of 7 is equal to 
16, which is not significantly different from the value of 18 in equa- 
tion (17). 

Dyer [1] has analyzed experimental estimates of 0 q in the Monin 
layer and finds that 


0 Q = (1 - 15z/L 0 )-°‘ 55 . 


( 21 ) 


This result is not significantly different from equation (17). Prasad 
and Panofsky [3] find that equation (17) summarizes the temperature obser- 
vations in their analysis. However, the largest value of -z/L 0 in these 
studies was on the order of one. It can be seen from the results in Fig- 
ure 1 that it would be rather difficult to find a dependence of 7 on R L 
if we restricted the analysis to the data points associated with 
-z/L q <1.0. The large values of -z/L Q in this study were primarily 
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obtained by making temperature measurements at sufficiently great heights 
(150 meter tower). Dyer [1] obtained his temperature measurements at 
various height intervals up to a maximum height of 16 m. The tempera- 
ture data used by Prasad and Panofsky [3] were obtained from various 
tower sites at heights between the 1- and 46-meter levels. In view of 
the result given by equation ( 20 ), we tentatively conclude that the 
result given by equation (17) for the Monin layer is not an asymptotic 
result as Rl approaches zero, but rather, equation (17) yields values of 
00 associated with typical values of — 0.004, 

Figure 8 is a plot of the experimental values of 0 q as a function of 
7z/L q , where y is calculated with equation (20). The function given by 
(19) is a good fit to the data, and a comparison of Figures 1 and 8 shows 
that the experimental scatter is reduced significantly by permitting the 
quantity 7 to be a function of R^. 

Upon combining equations (19) and (20), we find that 

0 e (z/L o ,R L ) - (1 - 0.0044R£ s/2 z/L 0 )‘ l/2 . (22) 

A plot of 00 as a function of z/L 0 for various values of R^, according to 
equation (22), is shown in Figure 5. For sufficiently large -z/L 0 or 
small R^, we have the asymptotic behavior 


0q(z/L o ,R l ) ~ (0.0044) -1 / 2 R3/ 4 (-z/L o ) -1 / 2 . (23) 


A comparison of Figures 1 and 5 will show that the expression given by 
( 22 ) is a reasonably good fit to the data. 

Equation (22) can be cast into two alternate forms, namely, 

0 e (z/L o ,fz/u* o ) = (1 + 0.0044(fz/u* o )" 3 / 2 (-z/L o ) 5 / 2 ) -1 / 2 

* 

(24) 

and 

0 Q (fz/i W R L ) = (1 + 0.0044 R £ s/2 fz/u * 0 )" l/2 . (25) 
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These functions are given in Figures 6 and 7. For sufficiently large -z/L Q 
or small fz/u, vo , we obtain from (24) the result 


0 0 ( Z /L o ,fz/u* o ) ~ (0.0044 r 1 / 2 (fz/u* 0 )3/ 4 (-z/L 0 )- 5 / 4 , (26) 


while for sufficiently large fz/u y , Q or small R^, we obtain from (25) 


0q (f z/ u^ 0 , R^) - (0.0044)- 1 / 2 R^/ 4 (fz/u^)- 1 / 2 . (27) 


These results are consistent with the experimental results in Figures 2 
and 3. 


V. POTENTIAL TEMPERATURE PROFILE 


Integration of the differential equation (19), subject to the boundary 
condition (5), will yield the potential temperature profile for any partic- 
ular value of R^, so that 


A9(z/z ) r( 1 + Az/ z ) 1 / 2 - 1 (1 + A) 1 / 2 + l~i 

In J V , (28) 

*° Hi + Az/z 0 ) 1 / 2 + 1 (1 + A) 1 / 2 - 1^ 

where 

A =- 7 z 0 /L o (29) 

A0(z/z o ) = 0(z/z o ) - 0(1) . (30) 


Figure 9 depicts A^5(z/z 0 )/T^ 0 as a function of z/ z Q for various values of 
A. As oo, the quantity A 0 and the function (29) has the asymptotic 

behavior 


A0(z/z o ) 


Svo 


In j- - 0.0011R£ 5/2 


Ro'Mf- - 1), 


(31) 
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where Ro is the drag Rossby number given by 


Ro = 


;-o /fz cr 


( 32 ) 


Thus, for sufficiently large R^, we have a logarithmic potential tempera- 
ture profile with a linear defect. 

An estimate of the potential temperature drop 60 across the planetary 
boundary can be obtained by evaluating equation (28) at the top of the 
boundary layer z = u Vc /4f, so that 


66 


= In 


*o 


(1 + 7/41^) x / 2 - 1 (1 - 7 z o /L q ) V 2 + 1 

(1 + 7 / 4 R ^) 1 / 2 +1 (1 - 7 Z 0 / L q ) 1 / 2 - 1 


(33) 


We can also express this result in the form 


66 


= In 


*0 


{ 


(1 4- 7/4R ^)^/ 2 - 1 (14- 7 /R^Ro ) 1 / 2 4- 1 

. 1, — n < — .■..i. w i.i.Mi m ■■■ 1-1.1 . » " « " »■ 1 . mmmmmmm r mm 

(1 + 7 / 4 RL ) 1 / 2 + 1 (1 + 7 / RjRo ) 1 / 2 - 1 


(34) 


As R l approaches zero, the quantity 7 approaches zero, so 
that upon expanding (34) in series, we find that 


66 

T 

Svo 


~ In - 0.0022 RT 5 / 3 (7 - Ro' 1 )+ ... . 


(35) 


This result shows that for sufficiently near neutral atmospheres, or 
atmospheres with sufficiently large rotation rates (f), the dimension- 
less potential temperature drop across the planetary layer is asymp- 
totically equal to ln(Ro/4) as R^ approaches infinity. 

As z 0 /L 0 becomes small, equation (33) behaves like 



- lD< 4 L? 


+ In 


(1 + 7 /4R l ) 1 / 2 - 1 

(1 + 7/4R l ) 1 / 2 + 1 


(36) 
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For sufficiently small R^, the second term on the right-hand side of 
(36)> is negligibly small, so that 

I 2 - - ln(227.3R 5 /2Ro). (37) 

This corresponds to the case in which the boundary layer is very unstable. 
If we would have neglected Coriolis effects ab initio and set 7 = 18, then 
the result that corresponds to (37) is 

f 2 - ~ -ln(-4.5z /L Q ), (38) 

vcO 

or in terms of and Ro, 

I 2 - - ln(0.222R L Ro). (39) 


VI. SECOND- AND THIRD- ORDER CONSIDERATIONS 

To facilitate the discussion that follows, let us denote 0 U and 0q 

as given by (.16) and (17) with 0 ^ and 0^^ and the associated value 

(l ) 

of 7 (=18) with 7 ^ . The calculations of the scaling temperature T 

and scaling velocity were based on (see Appendix A). However, 

^ (2) 

the experimental estimates of 0^ in Figure 1, 0^ say, implied that 7 

is not a constant equal to 18, but rather that 7 is a function of R~ , and 

(2) L 

we will denote this function with 7 ' . Th.us, in a manner of speaking, 

we have deduced the second-order function 0^ with the aid of the first- 

(1) ® 

order Monin layer through j2r ' . According to Businger’s hypothesis, 

0 ~ in the Monin layer; therefore, we might infer that in the 

Monin layer the second-order dependence of 0 U on resulting from the 
second-order potential temperature profile is 

0^ 2) = (1 - 7 (2) z/L 0 ) -1 / 4 . (40) 
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This result will be used to obtain new estimates of T* 0 and u* 0 , and 
thus third-order estimates of 0q and 7. Henceforth, a numerical 
parenthetical superscript will denote the order of the quantity in 
question. 

To analyze the effects of the third-order corrections on 7, it 

is useful to write in the form 

'u 

0^ = (1 - 18z/l/ 2) ) -l/4 , (41) 

where 

L^ s) = 18L o / 7 (2) . (42) 


The wind profile that corresponds to (41) can be written in the form 

„<2) 


Q (2 ) (z) = jin f- - 1 (r (2) (z/l/‘ 2) ) , 

it o 


(43) 


where 


-z/L 


(2) 


♦ ( 2 ) (z/L< 2 >) 


-/ 


1 - ^ 2 ) (i) 




(44) 


and | is a variable of integration associated with z/L^ ' . The ratios 
between the first- and second-order scaling parameters are 


“io* _ ln(z/z Q ) - ^ (l) (z/L^ i:> ) 
ln(z/z ) - i|/ 2 ^ (z/L^ 2 b 


(45) 


c (2) 

Svo 

tU) 

* o 


ln(z/ z Q ) - i|/ 1 ^(z/l/ 1 b'| 2 
•ln(z/z 0 ) - t (2) (2/L,$ 2) ) 


(46) 
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where (z/Lq 1 ^ ) is the first-order counterpart of i|/ 2 ^ (z/l/ 2 ^ ) * The 

/ \ / \ o 

18-meter level wind speed was used to obtain u^ Q and ^*0 > so that 

( 2) 

z = 18 m and z Q 2 0.2 m at the NASA 150-meter tower site. If 7 = 18, 
then z/l/ 2 ^ - z/L q , and it follows from (45) and (46) that u^ 2 ^ = u*q^ 
and T^ 2 ^ - T^\ ^ 7 ^ is ^ ess t ^ ian 18, then -z/l/ 2 ^ is less than 

-(z/L^ 1 ^), which implies -\j/ 2 ^ (z/L^ 2 ^ ) < (z/L^ 1 ^ ) , and we conclude 

from (45) and (46) that u^ 2 ^ is less than and |T^ 2 ^ | is less than 

|T*o^l* ^ I s greater than 18, the inequalities are reversed so 

that u£ 2 ^ is greater than u^ 2 ^ and | | is greater than |T^ 2 ^|. At 
z/L = 0 , we have uf^ = ufit and and the ratios 

and de Part from unity as -z/L 0 approaches infinity. Now 0 is 

(2) (3) (2) 

directly proportional to T ^ Q 7 and 7 is inversely proportional to u J j Q / . 

(3) 

This means that the third-order estimates of 7 , 7 , will lie below and 

(2) (2) 

to the right of their corresponding estimates of j J for y K 7 < 18 and 

( 2) 

above and to the left of their corresponding estimates, of 7 for 
(2) 

7 > 18. The third-order estimates of 7 are shown in Figure 10. A 

least-squares fit to the data points in this figure yields the result 

7 (3) = 0.0012(R^ 2 V 7 / 4 . 

This function is depicted in Figure 9 along with the functions 7 ^ and 
( 2 ) 

7 V 7 . The results in the figure show that this correction procedure seems 
to converge. The differences between 7 ^^ and 7 ^ are probably not 
significant because 7 ^ lies well within the data scatter of 7 ^ and 
vice versa. 

The question that must now be answered is whether or not the func- 
(2) 

tion j 2 r as expressed by (40) is physically real or just a spurious 
result. This can be answered only by analyzing wind and temperature 
profile measurements with corresponding direct measurements of the sur- 
face tangential stresses and heat fluxes. 
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VII. CONCLUDING COMMENTS 


We have analyzed three representations of 0 0 in the planetary 
boundary layer, namely, 0 0 (z/L o ,Rl), 0q(z/L o , fz/u* 0 ), and 0 0 (fz/u* o , 

R L ) . However, it is reasonable to expect that 0 0 should depend on 
other parameters. For example, it is well known that baroclinic 
effects manifested by the pressure gradient terms in the horizontal 
momentum conservation equations produce significant effects on the 
wind profile [11], and it is not unreasonable to assume that these 
effects are also reflected in the potential temperature profile. The 
available results of investigations of the wind profile in the planetary 
boundary layer could be used to suggest those additional parameters that 
should be included in future attempts to model the potential temperature 
profile. 

In our analysis we developed an empirical model of 0 0 (z/L o ,Rl), and 
from this function we deduced the functions 0 0 (z/L o , fz/u^ Q ) and 
0 q ( fz/u^QjR^) which agreed reasonably well with the experimental results 
However, in studies of this kind, one must be careful to avoid or to at 
least recognize the problem of the spurious correlation which results 
when a scaling parameter occurs in two* or more dimensionless quantities. 
In this report, the scaling temperature T* 0 was derived from the esti- 
mates of L q and u. Vo , so that 0 0 is proportional to zL Q u£§ Sg/Sz. Thus, 
by the very fact that (1) z and L Q are contained in 0 0 and z/L 0 , (2) u* 0 
and L 0 are contained in 0 0 and R^, and (3) z and u* Q are contained in 0 0 
and fz/u. Vo , the dependencies of 0 0 on (z/L 0 ,R;l)“, (z/L q , f z/u* 0 )- , and 
(fz/u* 0 ,R;L)-coordinates implied in Figures 1 through 3 could be spurious 
This problem can be resolved by increasing the sample size and applying 
appropriate statistical significance tests. 

To verify the results in this report, simultaneous observations of 
the potential temperature profile, the surface heat flux H 0 , and the 
surface tangential stress r Q are required. The measurements of the 
potential temperature profile should extend up to sufficiently great 
heights in the Ekman layer (z ~ 0 (u^ Vo /4f ) ) . * The problem of spurious cor- 
relation between 0 0 and the various independent variables will still 
remain, because T* 0 is proportional to H q t " 1 / 2 and L 0 is proportional 
to so that 0 0 , z/L 0 , R l and fz/u^ Q are proportional to 


zT 1 / 2 H" 1 Se/Sz, zt~ 3 / 2 H , t H* 1 , and zt” 1 / 2 , 
o o 5 oo 5 o o 5 o 9 

respectively. The quantity t q is contained in all four of these quanti- 
ties; H Q is present in 0 0 , z/L 0 , and R L ; and the height z is contained 
in 00, z/L 0 , and fz/.u* Q . 
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With these comments in mind, we 
potential temperature profile in the 
influenced by the action of Coriolis 


provisionally conclude that the 
planetary boundary layer is 
forces on the mean flow* 
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Figure 1. Experimental Values of 0 , 
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Figure 2. Experimental Values of as a Function of z/L Q and fz/u* Q 
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Figure 3* Experimental Values of 0q as a Function of fz/u^ Q and 
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Figure 5. The Quantity 0 q as a Function of z/L 0 and 
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Figure &. The Quantity 0 q as 
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APPENDIX A 


Calculation of the Scaling Parameters u, Vo and T*o 


In the unstable Monin layer, the dimensionless mean flow shear is 
a universal function of z/L Q , so that 


k i z Bu _ . 

U, vo dz 4< z/ 


(A-l) 


where u(z) is the mean wind speed at height 2 , k x is von Karman’s con- 
stant with numerical value approximately equal 0.4, u. vo is the surface 
friction velocity, and 0 u (z/L o ) is a universal function of z/L 0 . The 
quantity L 0 is the surface Monin-Obukhov stability length, namely, 


L 


o 


4o C pP T 

k 1 g Ho 


(A- 2) 


In this equation H 0 is the surface heat flux, p and T denote the mean 
flow density and Kelvin temperature, g is the acceleration of gravity, 
and Cp is the specific heat at constant pressure. The dimensionless 
shear is related to the flux Richardson number through the experi- 
mentally derived relationship 


0 u = (1 - 18RiW\ 


(A- 3) 


which is given in reference 6. The flux Richardson number is defined as 

s. Se 

0 Bz 

R . = , (A- 4) 

Qu/bz ) 2 


where 0 is the mean potential temperature at height z. The flux 
Richardson number is a function of z/L Q . We shall invoke the Businger 
hypothesis [10] to relate Ri to z/L OJ so that 


Ri = z/L 0 . 


(A- 5) 
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Upon combining equations (A-l), (A-3), and (A-5) and integrating the 
resulting relationship, we find that 


u(z) 



where 


(A- 6) 





i - a + mnh dg 


(A- 7) 


We have used the condition that u(z Q ) = 0 in the derivation of (A-6), 
where z Q is the surface roughness length. Equation (A-7) can be eval- 
uated numerically for any value of z/L 0 . The lower bound of this 
integral may be set equal to zero because the contribution to from 
the region 0 < ~z/L 0 < -z 0 /L 0 is negligibly small. The surface fric- 
tion temperature T* 0 is defined as 


T 


*o 


1 

k l u *o 



(A- 8) 


Equations (A-2) and (A-4) through (A-8) can be used to calculate the 
scaling velocity u* 0 , the scaling temperature T Vc0 , and the surface Monin- 
Obukhov stability length. The procedure for calculating these quantities 
is as follows: (1) calculate the gradient Richardson number, equation 

(A-4), with the mean flow wind and potential temperature profile data in 
the Monin layer, (2) calculate L 0 with equation (A-5), (3) calculate u* Q 
with (A-6) and (A-7), and (4) finally calculate T*o with the results of 
steps (2) and (3) and (A-2) and (A-8). 

The 18- and 30-meter temperature and wind data were used to esti- 
mate the Richardson number. The distributions of u(z) and 0(z) between 
the 18- and 30-meter levels were assumed to be logarithmic profiles. The 
expressions for u(z) and 0(z) were differentiated with respect to z and 
evaluated at the geometric height z = 23 meters to yield estimates of 
3u/c)z and do/dz for the calculation of Ri. This Richardson number was 
used to calculate L 0 in step (2) above. The 18-meter level wind speed 
and \|r(l8m/L 0 ) were used to calculate u* Q in step (3). 

The surface roughness lengths z 0 that were associated with the NASA 
150-meter meteorological tower site and that were used in the calculation 
of u* 0 are given in reference 9. 
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APPENDIX B 


Wind Speed and Temperature Profile Data and Other Parameters 


The wind profile and. temperature data that were used in this report 
are given in Tables B-l and B-2. The 18-meter level wind direction is 
tabulated here because the surface roughness length at the NASA 150-meter 
tower facility is a function of wind direction. The values of Ri(23 m) , 
L OJ u *o> an< ^ T*o that were calculated with these data are given in 
Table B-3. 


TABLE B-l 


Table of Temperature Profile Data* 
(Temperature in °F) 


Case No. 

Date 

Time 

(EST) 

T(3m) 

(AT) x 

(AT) 2 

(AT) 3 

(AT) 4 

(AT) 5 

299 

1/23/68 

1315-1400 

76.0 

-1-.41 

-1.73 

-2.56 

-3.94 

-4.58 

305 

1/26/68 

1130-1230 

52.4 

-2.19 

-2.66 

-3.59 

-5.01 

-5.68 

310 

2/8/68 

915-1015 

40.5 

- .32 

-1.82 

-2.86 

-4.45 

-5.09 

319 

2/26/68 

1030-1130 

59.0 

-1.28 

-2.63 

-3.56 

-5.01 

-5.69 

323 

2/27/68 

1200-1300 

62.0 

-1.80 

-2.61, 

-3.57 

-5.36 

-5.95 

332 

3/1/68 

1210-1310 

49.0 

-2.31 

-2.94 

-4.04 

-5.73 

-6.54 

355 

3/20/68 

1330-1430 

75.0 

-1.82 

-2.21 

-2.99 

-4.32 

-4.99 

359 

3/22/68 

1200-1300 

77.0 

-1.92 

-2.97 

-4.33 

-5.57 

-6.18 

361 

3/22/68 

1320-1332 

78.0 

-2.00 

-2.83 

-4.17 

-5.65 

-6.01 

364 

3/27/68 

1408-1438 

73.0 

-1.89 

-2.32 

-3.39 

-4.86 

-5.53 

365 

3/28/68 

905-932 

73.5 

-1.44 

-2.42 

-3.12 

-4.52 

-5.26 

366 

3/28/68 

1145-1245 

74.5 

-2.63 

-3.13 

-4.09 

-5.62 

-6.37 

367 

3/29/68 

1010-1040 

73.8 

-1.45 

-1.96 

-2.61 

-4.20 

-5.18 

406 

4/7/68 

1507-1537 

79.0 

-2.02 

-2.57 

-3.55 

-4.88 

-5.48 

445 

5/7/68 

1400-1415 

78.2 

-2.29 

-2.54 

-3.74 

-5.33 

-5.99 

515 

6/16/68 

1200-1300 

83.0 

-2.16 

-2.60 

-3.65 

-5.31 

-5.98 

551 

6/29/68 

947-1027 

79.8 

- .93 

-1.90 

-2.48 

-3.95 

-4.52 

554 

6/30/68 

933-952 

83.1 

-1.40 

-1.91 

-2.67 

-4.04 

-4.67 

625 

10/14/68 

1100-1200 

80.0 

- .92 

-1.91 

-2.82 

-4.51 

-5.28 

i 


*(£0! = T(18m) - T(3m) 
(AT) 2 = T(30m) - T(3m) 
(AT) 3 - T(60m) - T(3m) 
(AT) 4 = T(120m) - T(3m) 
(AT) 5 = T (150m) - T(3m) 
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APPENDIX C 


Correction of Finite Difference Estimates 
of Vertical Gradients of Potential Temperature 


A first estimate of Sg/Sz can be obtained with finite differences, 
so that 


oe/5z) (1) 


§2 " §1 
z 2 - Zi 5 


(C-l) 


where the superscript on Sg/Sz denotes the order of the estimate, and 
0 X and Q 2 are the potential temperatures at heights z 1 and z 2 . A 
second and more precise estimate of <$Q/dz can be obtained by assuming 
that Sq/Bz is related to z through a power law 


dQ/dz = az q , 


(C-2) 


where a and q are constants. This assumption is appropriate if it is 
applied in a piecewise manner. The quantity q will vary between -1 and 
-2. Upon integrating (C-2) between levels z x and z 2 , we find that 


§2 


01 ” 


a 

i + q 



(C-3) 


Elimination of h between (C-2) and (C-3) and evaluation of the result- 
ing relationship at the midpoint z = (z x + z 2 )/2 of the interval 
z x ^ z ^ z 2 yield a second estimate of S§/Sz, so that 


/ \ (1 + q) (02 " 0i) ( z 2 z i) 

(Se/32) ' 2’ ^ ‘ 


(C-4) 


Upon forming the ratio between (C-l) and (C-4), we find 


r 


0, 


( 2 ) 


0, 


(i) 


(1 + q) Cl - l)(l + l) q 
2 q (l - | 1+q ) 


(C-5) 
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where 



(C-6) 

(C-7) 

(C-8) 


Equations (C-l), (C-5), and (C-7) were used to estimate the dimension- 
less temperature gradient 0q. The procedure consisted of the follow- 
ing steps: 

(1) Obtain a first estimate of d0/ciz with (C-l). 

(2) Calculate a first estimate of the dimensionless potential 
temperature gradient 0Q 1 ) at z = (z x + z 2 )/2 with (C-7), 

(3) Plot 0Q 1 ) as a function of -z/L 0 on bilogarithmic graph 
paper and estimate piecewise values of q. 

(4) Calculate the quantity r with (C-5) using the results of 

(3). 

(5) Calculate a second estimate of 0 q by multiplying the 
yalues of 0 q 1 ^ obtained in step (2) with r. 

The results of step (5) above were used to model the potential tempera- 
ture profile. 

A plot of r as a function of q for | = 0.167, 0.5, 0.6, and 0.8 is 
shown in Figure C-l. These values of r correspond to the ratios between 
heights appropriate for the NASA 150-meter meteorological tower site. 
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